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ABSTRACT 

We determine the orbits of four double degenerate systems (DDs), composed of two 
white dwarfs, and of two white dwarf - M dwarf binaries. The four DDs, WD1022+050, 
WD1428+373, WD1824+040, and WD2032+188, show orbital periods of 1.157155(5) 
d, 1.15674(2) d, 6.26602(6) d and 5.0846(3) d respectively. These periods combined 
with estimates for the masses of the brighter component, based on their effective tem- 
peratures, allow us to constrain the masses of the unseen companions. We estimate that 
the upper limit for the contribution of the unseen companions to the total luminosity in 
the four DDs ranges between 10 and 20 per cent. In the case of the two white dwarf - M 
dwarf binaries, WD1042— 690 and WD2009+622, we calculate the orbital parameters 
by fitting simultaneously the absorption line from the white dwarf and the emission 
core from the M-dwarf. Their orbital periods are 0.337083(1) d and 0.741226(2) d 
respectively. We find signatures of irradiation on the inner face of WD2009+622's 
companion. We calculate the masses of both components from the gravitational red- 
shift and the mass-radius relationship for white dwarfs and find masses of 0.75 - 0.78 
M Q and 0.61 - 0.64 M for WD1042-690 and WD2009+622 respectively. This indi- 
cates that the stars probably reached the asymptotic giant branch in their evolution 
before entering a common envelope phase. These two white dwarf - M dwarf binaries 
will become cataclysmic variables, although not within a Hubble time, with orbital 
periods below the period gap. 
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1 INTRODUCTION 

About 10 per cent of white dwarfs reside in close binary 
systems (Napiwotzki et al. 2003). The formation of white- 
dwarf close-binary systems commences with a pair of main- 

* The Isaac Newton and William Herschel telescopes are oper- 
ated on the island of La Palma by the Isaac Newton Group in the 
Spanish Observatorio del Roque de los Muchachos of the Instituto 
de Astrofi'sica de Canarias. Based on observations collected at the 
Centro Astronomico Hispano Aleman (CAHA) at Calar Alto, op- 
erated jointly by the Max-Planck Institut fur Astronomie and the 
Instituto de Astrofi'sica de Andaluci'a (CSIC). Based on observa- 
tions collected with ESO telescopes at the Paranal Observatory 
under programme IDs 165.H.-0588 and 167.D-0407. 



sequence stars orbiting in a wide binary. The most mas- 
sive star will evolve faster becoming a giant and transferring 
mass to its companion. If mass transfer is sufficiently rapid, 
the result will be the formation of a common envelope (CE) 
made of the outer layers of the giant. The orbital energy of 
the binary can then be used to eject the envelope, the re- 
sult being a binary composed of a white dwarf and a main 
sequence companion. If the mass ratio of the initial compo- 
nents was close to unity, the binary formed after CE ejection 
can still be wide, whereas in the case of extreme mass ra- 
tios, the two components will spiral in to eject the envelope 
forming a tighter binary. This binary, called a post common 
envelope binary (PCEB), if close enough, might become a 
cataclysmic variable (CV) if mass is transferred stably from 
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the main sequence star to the white dwarf. If, on the other 
hand, the binary is too wide to become a CV, the main 
sequence star will evolve into a giant and the system will 
undergo another CE phase. When the envelope is ejected, 
the resulting binary will be composed of two white dwarfs a 
few solar radii apart - called a double degenerate (DD). For 
details on the evolution and formation of DDs see Nelemans 
et al.'s (2001) recent population synthesis studies and ref- 
erences therein. The study of white-dwarf close-binary stars 
can help us understand the elusive CE phase that they must 
have gone through, at least once, during their evolution. This 
phase is very difficult to study in any other manner as it lasts 
only of the order of 1 to 100 years. For recent calculations 
on CE ejection efficiencies see Soker & Harpaz (2003) and 
references therein. 

The subjects of this paper, white dwarf - M dwarf bi- 
naries and DDs, are, as mentioned above, the progenitors 
of CVs, and potential progenitors of Type la supernova re- 
spectively. A short period DD where the combined mass of 
the white dwarfs exceeds the Chandrasekhar mass might 
become a Type la supernova (Iben & Tutukov 1984) , but 
no candidates have yet been found. We should mention that 
this idea is still quite controversial with some theorists claim- 
ing that a DD does not become a Type la supernova (Saio 
& Nomoto 1998) and others claiming that this can be the 
case if one includes rotation in the calculations (Piersanti et 
al. 2003). The ESO Supernova la Progenitor Survey (SPY) 
has as its main goal to search systematically for massive, 
short period DDs in the galaxy to establish their direct link 
with Type la supernova and it is the source of more than 100 
recently discovered DD systems (Napiwotzki et al. 2003). Ef- 
forts are also been carried out to systematically find white 
- M dwarf binaries by using the Sloan Digital Sky Survey 
(SDSS) resulting in more than 400 possible candidates (Sil- 
vestri, Hawley & Szkody2003). 

Schreiber & Gansicke (2003) review the evolution of 
the known sample of white-dwarf close binaries with main 
sequence companions, the PCEBs, and conclude that, due 
to selection effects, the (until then) known population of 
PCEBs is dominated by young systems with hot white 
dwarfs that will evolve into short period CVs (P < 3 h). 
The PCEB sample is biased toward low mass companions 
as well which is the reason why they evolve into short period 
CVs. In contrast, the PCEBs found in the SDSS (Raymond 
et al. 2003), published later on, show an average white dwarf 
temperature significantly lower, demonstrating that to sam- 
ple the full parameter space better selection criteria have to 
be devised. 

In this paper we study 6 white dwarf close-binary stars. 
Four of them are found to be DDs, the other two are com- 
posed of a white dwarf and an M dwarf star. The separation 
of both components in these binaries is of the order of a 
few solar radii, as determined from their short orbital peri- 
ods, so they are detached systems where no mass transfer 
between them takes place. We obtain their orbital solutions 
and compare the results obtained with predictions drawn 
from population synthesis studies. 



Table 1. Journal of observations. A description of the setup used 
with each telescope is given in the text. N indicates the number 
of spectra taken with each setup. * indicates high resolution spec- 
tra covering only the emission core of Ha. The brightness in B 
magnitudes of each target is also given. 



Object 


N 


Date 


Setup 


Inst. 


WD1022+050 


10 


29/2-3/3/96 


AAT 


RGO 


(B = 14.37) 


3 


19/3/97 


AAT 


RGO 




6 


8-10/2/98 


INTa 


IDS 




2 


5/6/98 


AAT 


RGO 




2 


15/4/03 


INTb 


IDS 




2 


17/4/03 


INTb 


IDS 


WD1428+373 


2 


16/2/97 


INTb 


IDS 




4 


y-J-U/ z I yo 


INTa 


IDS 




9 


3-5/3/99 


INTa 


IDS 




8 


23/7/00 


WHT 


ISIS 




2 


23/1/03 


WHT 


ISIS 


WD1824+040 


13 


19-24/6/95 


INTa 


IDS 


(B = 14.00) 


13 


29/2-3/3/96 


AAT 


RGO 




17 


1 8-91 /Q7 


AAT 


RGO 




2 


23-24/6/97 


INTa 


IDS 




2 


17-18/5/00 


VLT 


UVES 




5 


4/5-18/8/01 


VLT 


UVES 




1U 


zo-oU/ 1U/U1 


1IN la 


TTiC 

IJJo 




O 


oq 97 /o /no 


o. om 


1 VV UN 


WD2032+188 


10 


12-2A/Q/m 


WHT 


ISIS 


(B = 15.30) 


3 


14-15/8/93 


WHT 


ISIS 




8 


10-12/6/95 


WHT 


ISIS 




1 


6/11/97 


INTa 


IDS 




1 


22/11/97 


WHT 


ISIS 




2 


15/4/03 


INTb 


IDS 


WD1042-690 


15 


29/2-3/3/96 


AAT 


RGO 


(B = 13.05) 


5 


18-20/3/97 


AAT 


RGO 


WD2009+622 


4* 


10/6/96 


WHT 


UES 


(B = 15.00) 


6 


21-25/6/95 


INTa 


IDS 




6 


11/7/98 


WHT 


ISIS 




8 


14/7/98 


WHT 


ISIS 




7 


6-7/10/98 


WHT 


ISIS 



2 OBSERVATIONS AND REDUCTION 

The data used in this study were taken over many years 
(since 1993) using five different telescopes and seven differ- 
ent setups. Table 1 gives a list of the number of spectra taken 
for each target at each observing campaign indicating also 
which telescope was used in each case. 

AAT: denotes data taken with the Royal Greenwich Ob- 
servatory (RGO) spectrograph at the 4 m Anglo- Australian 
Telescope (AAT) . The setup consisted of the 82 cm camera 
with the R1200R grating centred in Ha. The CCD used was 
an MIT-LL (3kxlk) in fast readout mode. This combination 
gives a dispersion of 0.23 A pixel -1 . 

INTa: denotes data taken with the Intermediate Disper- 
sion Spectrograph (IDS) at the 2.5m Isaac Newton Telescope 
(INT) on the island of La Palma. For these data, the setup 
consisted of the 500 mm camera with the R1200R grating 
centred in Ha and the Tek (lkxlk) chip. This combination 
results in a dispersion of 0.39 A pixel -1 . 

INTb: denotes data taken also with the IDS at the 
INT but with a setup that consisted of the 235 mm camera 
with the R1200B grating and the EEV10 (2kx4k) CCD. The 



© 0000 RAS, MNRAS 000, 1-16 



Six detached white-dwarf close binaries 3 



wavelength range covered in this case included H7 and H/3. 
This combination results in a dispersion of 0.48 A pixel" 1 . 

WHT: denotes spectra taken with the 4.2m William 
Herschel Telescope (WHT) on La Palma. Most of the spec- 
tra were obtained using the double arm spectrograph ISIS. 
Only the red spectra were used in our study. Except for the 
spectra taken in January 2003, the setup for the red arm con- 
sisted of the 500 mm camera with the R1200R grating and 
a Tek CCD (lkxlk) giving a dispersion of 0.40 A pixel -1 . 
For the spectra of WD 1428+373 taken in January 2003, a 
MARCONI CCD (2kx4.7k) was used giving a dispersion of 
0.23 A pixel -1 . In the case of the 4 high resolution spec- 
tra taken of WD2009+622 (marked with a star in Table 1), 
the Utrecht Echelle spectrograph (UES) was used with the 
35 cross disperser and a Tek (lkxlk) CCD giving a typical 
dispersion of 0.07 A pixel -1 . 

VLT: denotes data taken with the UV- Visual Echelle 
Spectrograph (UVES) at the UT2 (Kueyen) VLT 8.2 m tele- 
scope located at Paranal Observatory. The setup used con- 
sisted of Dichroic 1 (central wavelengths 3900 A and 5640 A) 
with a EEV CCD (2kx4k) for the blue arm and two CCDs, 
a EEV (2kx4k) and a MIT-LL (2kx4k) for the red arm. This 
setup allows us to achieve almost complete spectral cover- 
age from 3200 A to 6650 A with only two ~80 A wide gaps 
at 4580 A and 5640 A. A slit width of 2.1" was used to min- 
imise slit losses and the CCDs were binned 2 x 2 to reduce 
readout noise. This setup results in a spectral resolution of 
0.36 A (or better if the seeing disk is smaller than the slit 
width) at Ha. Exposure times were 5 min. 

3.5 m: denotes data taken with the double beam TWIN 
spectrograph at the 3.5 m telescope in the Calar Alto Ob- 
servatory. Only the red spectra were used in this paper. 
The setup for the red arm consisted of the 230 mm camera 
with the T06 grating (1200 grooves mm -1 ) and a SITe-CCD 
(2kx0.8k) giving a dispersion of 0.55 A pixel -1 . 

The slit width was set to values between 0.8 and 1 arcsec 
depending on the seeing. We made sure in every case that 
the star filled the slit to avoid systematic errors in the radial 
velocities caused by the star wandering in the slit during an 
exposure. 

For the AAT, INT, and WHT runs we obtained CuAr 
plus CuNe frames to calibrate the spectra in wavelength. 
In the case of the 3.5 m and VLT observations the wave- 
length calibration arc used was a ThAr. All target spectra 
were bracketed by arc spectra taken within one hour and the 
wavelength scale interpolated to the time of mid-exposure. 
We subtracted from each image a constant bias level deter- 
mined from the mean value in its over-scan region. Tungsten 
Hatfield frames were obtained each night to correct for the 
pixel to pixel response variations of the chip. Sky flatfields 
were also obtained to correct for the pixel to pixel variations 
of the chip along the slit. After debiasing and flatfielding the 
frames, spectral extraction proceeded according to the op- 
timal algorithm of Marsh (1989). The arcs were extracted 
using the profile associated with their corresponding target 
to avoid systematic errors caused by the spectra being tilted. 
Uncertainties on every point were propagated through every 
stage of the data reduction. We did not attempt to correct 
for light losses in the slit. For the VLT spectra a special 
procedure was applied to correct for a quasi-periodic rip- 
ple pattern appearing in many of the uncorrected merged 
spectra (Napiwotzki et al. 2005) . The resulting VLT spectra 



were then divided by a smoothed spectrum of a DC white 
dwarf, which by definition shows no spectral features at all 
and therefore provides an excellent means of correcting for 
the instrumental response. 



3 RESULTS 

3.1 Average spectra 

Fig. 1 presents the average red spectra for the six sys- 
tems discussed in this paper. The spectra of four of the 
systems, WD1022+050, WD1428+373, WD1824+040 and 
WD2032+188, are very similar showing only very broad ab- 
sorption at Ha. WD1042-690 and WD2009+622, on the 
other hand, show a double line structure in their Ha line 
profile composed of broad absorption coming from the white 
dwarf and narrow emission from the heated surface of the 
M dwarf companion. As this narrow emission moves within 
the absorption profile with the orbital period, an average 
spectrum would show the emission broadened and for that 
reason for these two systems we present the spectrum at a 
particular orbital phase instead of an average of the spec- 
tra during an orbit. The narrow emission in WD 1042— 690 
is, on average, significantly stronger than for WD2009+622. 
For WD2009+622, the strength of the narrow emission core 
changes with orbital phase and only at orbital phase 0.5 
(when we are looking directly into the heated face of the 
M dwarf) reaches similar strength to that in WD1042— 690. 
See Section 3.3.1 for details. 



3.2 Four double degenerate systems 

To measure the radial velocities of the four double de- 
generates: WD1022+050, WD1428+373, WD1824+040 and 
WD2032+188, we used least squares fitting of a model line 
profile. The model line profile is the summation of three 
Gaussian profiles with different widths and depths. For any 
given star, the widths and depths of the Gaussians are op- 
timised and then held fixed while their velocity offsets from 
the rest wavelengths of the lines in question are fitted sep- 
arately for each spectrum; see Maxted, Marsh & Moran 
(2000c) for further details of this procedure. 

Once the radial velocities for each system were known 
(see Table 2) we used a "floating mean" periodogram to 
determine the periods of our targets (e.g. Cumming, Marcy 

6 Butler 1999). The method consists in fitting the data with 
a model composed of a sinusoid plus a constant of the form: 

7 + Ksin(27r/(t-t )), 

where / is the frequency and t is the observation time. The 
key point is that the systemic velocity is fitted at the same 
time as K and to. This corrects a failing of the well-known 
Lomb-Scargle (Lomb 1976; Scargle 1982) periodogram which 
starts by subtracting the mean of the data and then fits a 
plain sinusoid; this is not the best approach for small num- 
bers of points. We obtained the \ °f the fit as a function 
of / and then identified minima in this function. 

Table 3 gives a list of the orbital parameters derived for 
each DD binary star. The orbital period of the second best 
alias is also given, along with the difference in \ 2 between 
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Table 2. Radial velocities for the 5 systems. In the case of the 2 white dwarf-M dwarf systems the values given are measured from the 
M dwarf. 

HJD - 2440000 RV (km s" 1 ) HJD - 2440000 RV (km s" 1 ) HJD - 2440000 RV (km s" 1 ) 



WD1022+050 WD1824+040 WD2032+188 



10143 


.1257 


—A2.2A±2. 


.55 


9893.4920 


96.12± 3. 


76 


9162. 


5641 


42 


79±5 


58 


10143 


.1447 


o a no 1 o 

— 34.08±3. 


03 


9893.4999 


95.64± 4. 


97 


9162. 


5752 


52 


97±5 


,51 


10144 


.0893 


20.85±2. 


67 


10143.2756 


ro no 1 o 

58.98± 2. 


,83 


9162. 


6371 


54. 


49±4. 


,82 


10144 


.1012 


16.48±2. 


46 


10143.2818 


56.00± 3. 


,20 


9162. 


,6485 


53. 


30±5, 


,03 


10145 


.0577 


89.35±1. 


96 


10143.2880 


52.64± 2. 


,81 


9214. 


,5046 


97. 


,72±4, 


,18 


10145 


.0697 


81.96±2 


01 


10143.2949 


51.37± 3. 


,05 


9214. 


5220 


99 


70±4 


15 


10146 


.0863 


111.94±2. 


21 


"1 C\ "1 A A O 7^ C { • 

10144.2766 


"i nn a ^ 1 i 

100. 47± 1 


85 


9214. 


.5414 


102. 


,52±4, 


,47 


10146 


,0983 


116.99±2. 


04 


10144.2851 


100.91±1. 


96 


9878. 


.6788 


—7. 


75±4, 


,36 


10146 


.1719 


99.17±1. 


.99 


"1 C\ "1 A A O C\ O / ■ 

10144.2936 


102.73±2. 


10 


9878. 


,6848 


—5. 


,62±4 


46 


10146 


.1839 


1 \T A C\ I n 

97.42±2 


22 


10145.2760 


101.46±2 


65 


9879. 


7064 


62 


70±8. 


39 


10527 


.1487 


5.75±3. 


16 


10145.2892 


99.27±3. 


,30 


9879. 


.7136 


61 


,99±7 


84 


10527 


.1584 


8.34±3. 


63 


1 ri-l AC O Tn "1 

10146.2791 


58.95±2. 


.54 


9879. 


,7207 


72. 


,30±8 


,73 


10527 


.1670 


5.32±3. 


.83 


1 ri-i a c o o z' 1 n 

10146.2862 


f 77 1A 1 O 

57.30±2. 


88 


9880. 


7100 


98. 


97±7 


60 


10852 


.5812 


-31.10±4. 


68 


10146.2934 


57.87±3. 


59 


9880. 


7171 


95 


36±7. 


19 


10852 


.5911 


-41.51±3. 


72 


10146.2986 


54.33±8. 


,30 


9880. 


7270 


98 


67±7. 


,50 


10854 


.4604 


56.58±4. 


83 


10526.2345 


89.86±1. 


,74 


10759. 


3259 


57 


43±8. 


,35 


10854 


.4688 


66.99±4. 


17 


10526.2465 


92.02±1. 


.68 


10775. 


,3924 


97 


,21±3 


23 


10855 


.4728 


109.75±3. 


07 


10526.2585 


96.23 ±1 


81 


12745. 


,6895 


—39 


02±7. 


,58 


10855 


.4853 


109.26±2. 


80 


i n 7 n / ■ o i r" 

10526.2715 


95.24±1. 


64 


12745. 


,7036 


—22 


83±6, 


,93 


10969 


.8494 


102.51±3. 


08 


10526.2835 


94.06±1. 


,60 


WD1042- 


-690 






10969 


.8567 


107.65±3. 


.09 


10526.2955 


92.74±1. 


.36 


10143. 


1591 


76 


49±0. 


,39 


12744 


.4589 


-30.08±4. 


.44 


"i n 7 n / ■ nn'70 

10526.3078 


n T n 1 o 

97.26±2. 


.34 


10143. 


1687 


76 


13±0. 


,36 


12744 


.4731 


-40.03±4. 


63 


10527.2658 


"i nn a r\ \ a 

109.40±4, 


04 


10143. 


,9950 


—56. 


70±0. 


44 


12746 


.5091 


26.01±3. 


96 


10527.2777 


"I 1 O C f I o 

113.65±3. 


27 


10144. 


0000 


—57. 


84±0. 


44 


12746 


.5232 


23.86±4. 


21 


10527.2898 


109.12±2. 


87 


10144. 


,1363 


57. 


57±0. 


,54 


WD1428+373 




10527.2995 


"i nn r\f\ I a 

109.09±4, 


06 


10144. 


1401 


61 


33±0 


56 


10495 


.7128 


-40.29±7. 


71 


1 nrrio c\c c a 

10528.2664 


63.33±3. 


97 


10144. 


2244 


54 


00±0 


48 


10495 


.7230 


-29.57± 8. 


,36 


10528.2809 


69.52±3. 


61 


10144. 


,2586 


9 


,87±0, 


,54 


10853 


.7482 


-7.31± 3. 


,91 


1 n tr o n o o o o 


10.56±42 


73 


10144. 


,2625 


5 


78±0 


57 


10853 


.7658 


-10.22± 3. 


99 


"i n tr o n oni"i 


9.85±2. 


.60 


10144. 


2663 


-1 


43±0 


67 


10854 


.6999 


40.78± 3. 


,68 


"i n tr o n inner 

10529.3025 


5.22±3. 


,03 


10144. 


,2702 


-3 


85±0, 


,57 


10854 


.7709 


34.01± 3 


,39 


"i n 7 n c\ o n n 7 

10529.3095 


12.36±5. 


,84 


10145. 


,0378 


-58. 


52±0, 


,55 


11240 


.5750 


-79.23± 6. 


.07 




rn no 1 O 

50.08±3. 


68 


10145. 


.0416 


-56 


30±0 


47 


11240 


.6267 


-65.38± 3. 


68 


10623.7062 


—6.21±4. 


38 


10145. 


9757 


-29 


74±0 


58 


11240 


.7237 


-11.07± 3. 


80 


11681.7165 


35.50±1. 


.12 


10145. 


,9808 


-34 


51±0 


57 


11241 


.5861 


-81.31± 3. 


,88 


"1 "1 c o n o n o o 

11682.8933 


— 15.40±0. 


,79 


10526. 


,1651 


20. 


20±0, 


,47 


11241 


.6592 


-89.18± 3. 


.36 


i nn n o on j 

12033.8514 


— 13.50±1. 


,12 


10527. 


,1754 


21 


97±0, 


45 


11241 


.7579 


-54.58± 4. 


.47 


12078.7265 


6.4±0. 


63 


10528. 


1260 


77 


15±0 


,75 


11242 


.5826 


-57.44± 3. 


76 


12116.5747 


o i o n 1 n 

21.30±0. 


52 


10528. 


.1406 


69 


36±0. 


.92 


11242 


.7032 


-88.63± 3. 


89 


"l n 1 "I T 7 / \ o / \ 

12117.5939 


o "i nn 1 n 

81.90±0. 


56 


10528. 


,1558 


59 


51±1 


04 


11242 


.7749 


-92.46± 4. 


,03 


12139.5163 


"i - > on 1 n 

13.20±2. 


.44 


WD2009+622 






11749 


.3945 


-95.69± 1 


.98 


12209.3141 


— 12.86±1. 


,41 


9889. 


,6260 


-25. 


63±3. 


04 


11749 


,4086 


-88.57± 2. 


23 


12209.3241 


— 15.21±1. 


51 


9891. 


.6758 


50 
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Figure 1. Average spectra for the six systems studied in this paper. 



Table 2. Continued. 



HJD - 2440000 RV (km s" 1 ) HJD - 2440000 RV (km s" 1 ) HJD - 2440000 RV (km s" 1 ) 



WD2009+622 



11008.7131 
11008.7190 
11093.4111 
11093.4748 
11093.5206 
11093.5673 
11094.3644 
11094.4819 
11094.5644 



42. 
50. 
-76. 
-146. 
-194. 
-224. 
-232. 
-164. 
-50 



67±2.73 
81±2.59 
67±3.14 
25±5.21 
63±3.10 
30±3.29 
52±2.78 
87±2.82 
85±3.58 



Table 3. List of the orbital periods measured for the four double degenerate systems studied. To, the systemic velocity, 7, the radial 
velocity semi-amplitude, K, the reduced \ 2 achieved for the best alias, the 2nd best alias and the x 2 difference between the 1st and 2nd 
aliases are also presented. When calculating the \ 2 for both aliases we have added in quadrature a systematic error that results in a 
reduced \ 2 ~1 ( see text for details). The number of data points used to calculate the orbital period is given in the final column under N. 



Object 



HJD (T ) 
-2400000 



Period (d) 



7 (km/s) 



K (km/s) 



2nd best alias (d) A\ 2 



N 



WD1022+050 
WD1428+373 
WD1824+040 
WD2032+188 



51445.262(5) 
51579.64(1) 
51108.192(9) 
50947.07(5) 



1.157155(5) 
1.15674(2) 
6.26600(5) 
5.0846(3) 



39.05±1.19 
-21.46±1.62 
47.95±0.40 
35.11±1.52 



74.77±1.16 1.25 

67.90±1.68 1.41 

61.87±0.55 1.24 

63.50±1.59 0.66 



8.1580(4) 
1.22640(1) 
0.5449639(5) 
9.8267(6) 



40 
42 

3064 
45 



25 
25 
67 
25 
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Figure 2. Radial velocity curves for the four double degenerate systems. The data have been folded on the orbital period in each case. 
See Table 3 for the list of periods, radial velocity semiamplitudes and systemic velocities. Included in each panel is a plot of the residuals 
to the fit. The vertical scale on which the residuals have been plotted is twice the scale on which the radial velocities are plotted. 
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Figure 3. Each panel presents \ 2 versus cycles/day obtained after the period search was carried out. The frequency with the smallest 
X 2 corresponds to the orbital frequency of the system. For clarity we have also included an inset showing a blow up of the region where 
the best period is. The number of radial velocity measurements used for the period search calculations, n, is shown in each panel. 



the two best periods found. The large difference in \ 2 in- 
dicates that the second best aliases are not plausible. The 
resulting radial velocity curves (folded in the orbital period) 
are presented in Fig. 2 and the corresponding periodograms 
(x 2 versus orbital frequency) in Fig. 3. Each panel in the 
periodogram includes a blow up of the region in frequency 
where the minimum x 2 is found. 

In each case, we compute the level of systematic un- 
certainty that when added in quadrature to our raw error 
estimates gives a reduced x 2 ~ 1- By doing this we are 
considering the un-accounted sources of error such as true 
variability of the star or slit-filling errors that cause the poor 
fits of a few stars. Such errors are unlikely to be correlated 
with either the orbit or with the statistical errors we esti- 
mate, and therefore we add a fixed quantity in quadrature 
with our statistical errors as opposed to applying a simple 
multiplicative scaling to them. In all cases we use a mini- 
mum value of 2kms~ 1 corresponding to l/10 th of a pixel 



which we believe to be a fair estimate of the true limits of 
our data. The last column of Table 4 gives the value of sys- 
tematic uncertainty that we have added in quadrature in 
each case. 

We then calculate the probability of the true orbital 
periods being further than 1 and 10 per cent from the values 
we obtained - see Morales-Rueda et al. (2003a) and Marsh, 
Dhillon & Duck (1995) for an explanation of the method 
used to calculate these probabilities - and present them in 
Table 4. 

In all cases, the probabilities of the quoted periods being 
wrong are very low and we are certain that the values given 
in Table 3 correspond to the true orbital solution. 

In the cases where the probability of the orbital period 
being further than 1 and 10 per cent from our favoured value 
is the same, the significant probability lies within a very 
small range around the best period, with all the significant 
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Table 4. List of probabilities that the true orbital period of a 
system lies further than 1 and 10 per cent from our favoured 
value given in Table 3. Numbers quoted are the logarithms in 
base 10 of the probabilities. Column number 4 gives the value of 
the systematic uncertainty that has been added in quadrature to 
the raw error to give a reduced \ 2 ~ 1- 



Object 



1% 



10% 



systematic error 
(kins" 1 ) 



WD1022+050 
WD1428+373 
WD1824+040 
WD2032+188 



-7.66 
-9.59 
-1000 
-9.89 



-7.69 
-13.08 
-1000 
-9.90 



Table 5. The mass functions, f m , of the unseen components to- 
gether with the larger lower limits obtained by assuming i = 90° 
and by substituting in the mass function equation our determi- 
nation of the mass of the brighter component, Mi . 



Object 



/ ra (M Q ) Mi(M Q ) 



M 2 (M ) 
lower limit 



WD1022+050 0.050 0.389 0.283 

WD1428+373 0.038 0.348 0.233 

WD1824+040 0.154 0.428 0.515 

WD2032+188 0.135 0.406 0.469 



competition (i.e. next best alias) placed outside the 10 per 
cent region around the best alias. 



3.2.1 The unseen component of the binary 

By knowing the radial velocity semiamplitude of one of the 
components of the binary (the observable component), K, 
and the orbital period of the system, we can then calculate 
the mass function of the unseen component by using: 



fm — 



Ml sin 3 i _ PK j 
(Mi + M 2 ) 2 ~~ 2ttG ' 



(1) 



where the subscripts "1" and "2" refer to the brighter and 
the dimmer components respectively. The mass function is 
the lower limit of the mass of the unseen component. Ta- 
ble 5 gives the mass functions of the unseen components 
for the four DDs studied. In two cases (WD1824+040 and 
WD2032+188) the companion's mass function is greater 
than O.IMq which corresponds to the mass of a late M dwarf 
if it is a main sequence star. The masses of the brighter com- 
ponents of the systems have been measured by fitting their 
hydrogen line profiles to stellar atmosphere models using the 
tracks by Althaus & Benvenuto (1997) and can be substi- 
tuted, together with the assumption of the orbital inclina- 
tion of the system being 90° , in the mass function equation 
to give a larger lower limit for the masses of the unseen 
components. These revised lower limits (also presented in 
Table 5) are all greater than O.IMq which indicates that 
the unseen companions cannot be main sequence stars be- 
cause if they were we should be able to detect them (Marsh 
et al. 1995). The unseen companions must therefore be also 
compact objects, probably white dwarfs. 

We searched for the signature of the faint companions 
by shifting out the fitted radial velocity for each binary 



and then looking for differences in the line profiles at the 
quadrature phases (Marsh et al. 1995), i.e. 0.25 and 0.75. 
The spectra at quadrature phases were obtained by aver- 
aging the spectra contained in two separate phase ranges, 
i.e. the spectra in the range from 0.1 to 0.4 were averaged 
to obtain the phase 0.25 spectrum, and the spectra in the 
range from 0.6 to 0.9 to obtain the phase 0.75 spectrum. The 
results are plotted in Fig. 4. Any contribution from the com- 
panion white dwarf should be seen as an asymmetry in the 
line profile at phase 0.25 that is mirrored at phase 0.75 with 
respect to the rest wavelength (Marsh et al. 1995). None of 
the four systems show a clear asymmetry of this type in the 
line profiles. 

A second test that can be carried out to look for the 
faint companion consists in shifting out the fitted radial ve- 
locity off the spectra and creating a mean spectrum by com- 
bining all the shifted spectra, subtracting this mean spec- 
trum from the individual shifted ones and plotting the re- 
sulting spectra in a stack or trail. This aids the eye to iden- 
tify any leftover absorption moving with the binary orbit. A 
Doppler map (Marsh & Home 1988) can also be computed 
from these stack of spectra. Any orbital motion leftover in 
the spectra would appear in the maps as a small absorption 
region located in the V x — axis of the velocity map. We do 
not find any indication of the presence of the unseen com- 
ponent in either trails or Doppler maps in any of the four 
systems studied. 

3.2.2 How faint is the unseen companion? 

We explore the question of how small the contribution of 
the faint component of the system has to be so as not to 
be detected using the methods discussed in the previous 
section. 

To answer this question we create synthetic spectra that 
include the absorption corresponding to the brighter compo- 
nent of the system (in the form of three Gaussians), scaled 
to the measured value for each system, plus some extra ab- 
sorption moving opposite to it (also represented by three 
Gaussians) and determine for what percentage of brightness, 
relative to the bright component, we should be able to detect 
the faint component by looking at the spectra around the 
quadrature phases (Marsh et al. 1995). This method assumes 
that the companion stars have a spectrum similar to that of 
the brighter component. Although the fainter white dwarf is 
cooler and its Ha line will be less deep in its spectrum, this 
seems like a reasonable assumption as the companions are 
also white dwarfs, unless of course they are not DA white 
dwarfs. 

In the case of WD 1022+050 wc find that the bright- 
ness of the companion must be less than 10 per cent the 
brightness of the bright component for us not to detect it. 
The values we find for WD 1428+373, WD 1824+040, and 
WD2032+188 are 9 per cent, 23 per cent and 17 per cent 
respectively. These values translate into magnitudes for the 
companions, in the 6400-6700A region, that are respectively 
2.5, 2.6, 1.6 and 1.9 fainter than that of the bright com- 
ponents. Using the calculated absolute V magnitudes for 
the bright components (10.68, 10.35, 10.83 and 10.24, see 
the Discussion Section) and the broadband colour indices 
for pure hydrogen and log g = 8 stellar atmosphere models 
computed by Bergeron, Wesemael & Beauchamp (1995), we 



© 0000 RAS, MNRAS 000, 1-16 



8 L. Morales- Rueda et al. 



WD1022+050 



WD1428+373 




WD 1824+040 



WD2032+188 
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Figure 4. The spectra averaged around the quadrature phases for the four systems studied. In each case, the lower spectrum corresponds 
to quadrature phase 0.25 and the top one to phase 0.75. There is no clear asymmetry in the line profile at phase 0.25 that gets mirrored 
in phase 0.75 for any of the systems. This indicates that we cannot detect the faint companion of the systems. 



obtain upper limits for the absolute R magnitudes of the 
faint components of 13.3, 13.0, 12.5 and 12.2 respectively. 

3.3 Two white dwarf/M dwarf binaries 

In the case of WD1042-690 (aka BPM 6502) and 
WD2009+622 the spectra are composed of absorption lines 
that have their origin in the white dwarf plus emission cores 
that have their origin in the M dwarf. This extra emission 
component makes the measuring of the radial velocities more 
complicated, as the absorption coming from the white dwarf 
has its core filled by the M dwarf emission. A way to mea- 
sure simultaneously the radial velocities of both components 
is to use least squares fitting of a model line profile as in the 
previous case but this time using a model line profile that 
is the sum of four Gaussian profiles. Three of the Gaussians 
fit the absorption component and one fits the emission. The 
steps followed to carry out these complex fits consisted of: 1) 
fitting only the emission lines with a single Gaussian func- 
tion and obtaining the radial velocities associated to the 
emission line for each spectrum, 2) calculating the orbital 
solution for the emission lines by means of obtaining a peri- 
odogram from the radial velocities measured, 3) using this 



orbital solution to fix the orbit of the three Gaussians that 
will fit the absorption component of the lines and 4) obtain- 
ing the radial velocity semiamplitude and systemic velocity 
for the white dwarf by fitting all the spectra simultaneously 
with one emission and three absorption Gaussians. 

This method was easily applicable to WD2009+622 as 
the emission coming from the M dwarf is comparable to the 
Ha absorption core (see Fig. 1). In the case of WD1042— 690 
the emission component is very strong compared with the 
absorption core making the fitting of the data more difficult 
and the results obtained less accurate. For WD1042— 690 we 
find 7 very close aliases with very similar values of \ 2 ■ Ta- 
ble 6 gives a list of the orbital solutions for the 7 aliases. 
The solutions for K2, 72, Ki and 71 are consistent within 
the errors for all aliases. Previous studies of WD1042— 690 
(Kawka et al. 2000) result in an orbit solution in which the 
orbital period is consistent with our alias number 3 and the 
values for K2, 72, Ki and 71 are consistent with those from 
our 7 aliases. In Table 7 we present the orbital solutions for 
WD1042-690 and WD2009+622. The results presented for 
WD1042— 690 correspond to the first alias shown in Table 6. 
In Fig. 5 we plot the radial velocities measured for the emis- 
sion and absorption components together with the best fits 
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Figure 5. Orbital solution for WD1042— 690 and WD2009+622. Included in each panel is a plot of the residuals to the fit to the emission 
line component. The vertical scale on which the residuals have been plotted is larger than the scale on which the radial velocities are 
plotted. 



Table 6. Orbital solution for the 7 aliases found for 
WD1042-690. 



P(d) 


To 


72 




K 2 


X2 re 




-2450000 










0.337083(1) 


335.5553(7) 


7.93±0 


58 


69.31±0.69 


0.65 


0.337380(1) 


335.7239(7) 


7.83±0 


58 


69.24±0.69 


0.74 


0.336786(1) 


335.7239(7) 


8.03±0 


58 


69.38±0.69 


0.75 


0.337678(1) 


335.5550(7) 


7.73±0 


58 


69.16±0.69 


1.01 


0.336490(1) 


335.5557(7) 


8.14±0 


58 


69.44±0.70 


1.03 


0.337977(1) 


335.7238(7) 


7.64±0 


58 


69.07±0.68 


1.47 


0.336194(1) 


335.7240(7) 


8.26±0 


58 


69.49±0.70 


1.48 



Table 7. List of the orbital periods measured for the two white 
dwarf-M dwarf systems studied. To, the systemic velocity, 7, the 
radial velocity semi-amplitude, K, for both the white dwarf and 
the M dwarf, and the reduced x 2 achieved for the best alias are 
given The number of data points used to calculate the orbital 
period is also given. 

WD1042-690 WD2009+622 



N 
P(d) 
To (d) _ 
K WD (kms 1 ) 
1WD (kms -1 ) 
Km (kms -1 ) 
7m (kms -1 ) 

q=M M /M WD 



20 

0.337083(1) 
2450335.5553(7) 
-15.43±1.07 
44.08±0.75 
69.31±0.69 
7.92±0.58 
0.65 
0.223±0.018 



31/27 
0.741226(2) 
2450491.920(1) 
-43.81±2.11 
-64.38±1.80 
144.72±1.27 
-90.58±0.94 

1.12 
0.303±0.017 



given in Table 7. The error bars in the radial velocities are 
smaller than the size of the symbols used to plot the data. 
Notice that there are 4 extra points in the fit to the radial 
velocity of the M dwarf for WD2009+622. This accounts 
for the 4 high resolution spectra taken with UES that only 
covered the core of the Ha line. 



3. 3. 1 The M dwarf companions 

We binned the spectra for both binary systems into 10 phase 
bins and plotted them in Fig. 6 as a stack of spectra with 
orbital phase in the vertical axis. The advantage of present- 
ing the spectra in this way is that we can explore the line 
flux variations during an orbit. In the case of WD1042— 690 
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Figure 6. Trailed phase binned spectra for WD1042— 690 and WD2009+622. The orbit is plotted twice. The variable brightness in the 
emission line component, coming from the M dwarf companion, is clear in WD2009+622. 
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Figure 7. Flux modulation seen on the emission from the M 
dwarf component ofWD1042-690 (top) and WD2009+622 (bot- 
tom). 



the strength of the emission line does not vary significantly 
with orbital phase. This is not the case for WD2009+622 
where the line flux increases at orbital phases around 0.5. 
This phase corresponds to the white dwarf and the M dwarf 
being aligned with the line of sight, the white dwarf being 
closer to us. At this phase we are looking to the heated face 
of the M dwarf. 

To determine how the line flux varies with phase, we 
fitted the line profiles once more with four Gaussians, three 




2 4 6 1506 1506 15101634 1624.5 1625 1625.5 1626 
Orbital phase Orbital phase Orbital phase 

Figure 8. Flux modulation seen on the emission from the M 
dwarf component of WD1042-690 (top) and WD2009+622 (bot- 
tom). This time the data has not been folded in the orbital period. 



for the absorption component and one for the emission com- 
ponent but this time we allowed the height of the emission 
component to vary. The results are displayed in Fig. 7. The 
variation in flux is very significant in WD2009+622 and as 
we mentioned earlier it can be explained with emission com- 
ing from the irradiated face of the M dwarf, the side facing 
the white dwarf. This has important consequences for the 
calculation of the mass of the white dwarf carried out in 
the following sections as we have to account for distortions 
on the M dwarf to correct its radial velocity semiamplitudc. 
When the companion is heavily irradiated the emission line 



© 0000 RAS, MNRAS 000, 1-16 



Six detached white-dwarf close binaries 11 



used to measure its radial velocity will give us the value asso- 
ciated with the irradiated face, not the value for the centre of 
mass of the companion. Maxted et al. (1998) find that prob- 
ably as a result of optical depth effects, Balmer emission 
lines induced by irradiation are significantly broadened ren- 
dering measurements of radial velocities from fitting these 
lines more inaccurate. 

In the case of WD1042— 690, there is also some flux 
modulation present at a fainter level but peaking at phase 
0.6 rather than phase 0.5. This is probably associated to 
chromospheric activity of the M dwarf rather than irradia- 
tion. To make sure that this variability is not orbital we have 
also plotted the flux of the emission component versus un- 
folded orbital phase in Fig. 8. Wee see that for WD1042-690 
the variability observed is not phase dependent implying 
that is intrinsic to the M dwarf, probably chromospheric. In 
the case of WD2009+622 the variability observed is larger 
and it peaks at phase 0.5 confirming that it is due to irradi- 
ation of the inner face of the M dwarf by the white dwarf. 



Table 8. Summary of the parameters for WD1042— 690 and 
WD2009+622. BA indicates Benvenuto & Althaus (1998). Sec 
text for explanations on the different values given. 



Parameter 


white dwarf 


M dwarf 


WD1042-690 


fm(M ) 


0.0116(2) 


0.00013(1) 


M(M ) from Eq 2 


0.0174(9) 


0.0039(4) 


M(Mq) from q 


0.75(5)/0.78(5) 


0.1665(5)/0. 1735(5) 


M(M Q ) from BA 


0.75(7)/0.78(7) 




M(M Q ) if CO WD 


0.72 
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WD2009+622 


fm(M ) 


0.233(4) 


0.0065(5) 


M(Mq) from Eq 2 


0.40(3) 


0.12(2) 


M(M ) from q 


0.61(3)/0.64(3) 


0.1845(5)/0.1925(5) 


M(M ) from BA 


0.61(3)/0.64(3) 




M(M Q ) if CO WD 


0.59 






60/59 





3.3.2 The masses of the components 

Once the radial velocity semiamplitudes for both compo- 
nents have been measured, together with the orbital period 
we can use Eq. 1 to obtain the mass function for the white 
dwarf and the M dwarf in each case. If we combine Eq. 1 
with q = Mm/Mwd = Kwd/Km we obtain larger lower 
limits for the masses of both components: 



M W D = 



M M 



PKm(K wd + K m ) 2 
2-kG sin 3 i 
PKwd(Kwd + Km) 2 



(2) 



27rG sin 3 i 

The actual mass of the white dwarf can be determined 
from the gravitational redshift and the mass-radius relation- 
ship for white dwarfs (Althaus & Benvenuto 1997). First we 
calculate the difference in systemic velocities for both sys- 
tem components (jwd — Jm). We must then add corrections 
for (i) the redshift of the M dwarf, GMm / Rmc, where the 
radius of the M dwarf has been calculated by using the mass- 
radius relationship given by Caillault & Patterson (1990): 

logfl/R = 0.796 log M/M - 0.037, 

(ii) the difference in transverse Doppler shifts of both com- 
ponents: 

(Km - K 2 VD )/2csin 2 i, 

(iii) the potential at the M dwarf produced by the white 
dwarf: 

GMw d /ac, 

(iv) and the potential at the white dwarf due to the M dwarf: 

GM M /ac, 

where a is the distance between both stars. The mass of 
the white dwarf is then calculated by comparing the result- 
ing gravitational redshift with models for low mass helium 
white dwarfs (Althaus & Benvenuto 1997; Benvenuto & Al- 
thaus 1998). The inclination of the system can then be cal- 
culated by using Eq. 2. 

In order to calculate the corrections to the white dwarf's 



gravitational redshift we had to assume a value for the mass 
of the M dwarf. We performed several iterations of these 
calculations until we obtained consistent values for all the 
parameters. This method has been used previously (Marsh 
& Duck 1996; Maxted et al. 1998) to obtain the masses of 
both components in pre-CV systems. The results obtained 
after these iterations are presented in Table 8. 

We chose as the best estimates for the corrections de- 
scribed above and the final parameters of the iterations, 
those that resulted in a mass for the white dwarf that was 
closer to that calculated by Benvenuto & Althaus (1998). A 
minimum value for the masses of the components was found 
when we assumed a helium core white dwarf, with metallic- 
ity Z = 0.001 and an outer hydrogen envelope of fractional 
mass (i.e. mass of envelope/total mass of the star) 10 -s . We 
assumed a T cff = 21380 K and 25870 K for WD1042-690 
and WD2009+622 respectively (Bragaglia, Renzini & Berg- 
eron 1995; Bergeron, Saffer & Liebert. 1992). A maximum 
value was found if instead we assumed that the fractional 
mass of the outer hydrogen envelope was 4x 10 -4 in the case 
of WD1042-690 and 2xl0 -4 in the case of WD2009+622. 
The minimum and maximum values found for the masses 
are given in Table 8. 

For WD1042-690, jwD-fM is 36.16±2.72 kms -1 . Af- 
ter iterating for different values for the mass of the M dwarf 
we find that if Mm is 0.166 - 0.167 M the corrections i, ii, 
iii and iv are respectively 0.48, 0.09, 0.24 and —0.05 kms -1 , 
giving a value for the white dwarf redshift of 36.92 kms -1 . 
If, on the other hand, Mm is 0.173 - 0.174 M the correc- 
tions are respectively 0.48, 0.10, 0.25, —0.05 kms - , giving 
a value for the white dwarf redshift of 36.93 kms - . 

For WD2009+622, -/wd~Jm is 26.20±0.78 kms -1 . For 
a mass for the M M = 0.184 - 0.185 M or between 0.192 - 
0.193 M Q , the values for i, ii, iii, and iv are 0.49, 0.04, 0.12 
and —0.04 kms -1 , giving 26.82 kms -1 for the white dwarf 
redshift. 

If a carbon-oxygen core white dwarf with metallicity 
Z—0 and a hydrogen envelope of fractional mass 10 -4 is 
assumed instead, the values obtained for the masses of the 
white dwarf are 0.72 and 0.59 M for WD1042-690 and 
WD2009+622 respectively. 
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Table 9. New values for the M dwarf parameters depending on its filling factor / for WD2009+622. Rmcp is the radius calculated from 
the equation of Caillault & Patterson (1990) given in Section 3.3.2. A range of values for several parameters are given. These are the 
result of assuming two different masses for the outer hydrogen envelope of the white dwarf. See text for details. 



f 


K M 
km s _1 


q 


M M 

Mq 




i 

o 




a 




Rm 
R Q 




Rmcp 

Rq 




U.U 


1/1/1 70 
144. ( A 


U.oU 


n 1 c a 
U. lo4 


n i no 


an 
OU — 


oy 


o. iyi 


— O.Z04 


U.U 




U.zoy — 


U.Z4/ 


0.1 


149.95 


0.29 


0.178 


- 0.185 


63 - 


62 


3.182 


- 3.225 


0.085 - 


0.087 


0.233 - 


0.240 


0.2 


155.58 


0.28 


0.172 


- 0.179 


67 - 


65 


3.173 


- 3.216 


0.171 - 


0.173 


0.226 - 


0.233 


0.3 


161.64 


0.27 


0.165 


- 0.172 


72 - 


70 


3.165 


- 3.207 


0.256 - 


0.260 


0.219 - 


0.226 


0.4 


168.20 


0.26 


0.159 


- 0.165 


80 - 


76 


3.156 


- 3.198 


0.342 - 


0.346 


0.212 - 


0.219 


0.5 


175.30 


0.25 


0.152 


- 0.159 


79 - 


84 


3.147 


- 3.190 


0.427 - 


0.433 


0.205 - 


0.212 


0.6 


183.04 


0.24 


0.146 


- 0.152 


72 - 


74 


3.138 


- 3.180 


0.513 - 


0.520 


0.199 - 


0.205 


0.7 


191.49 


0.23 


0.140 


- 0.145 


67 - 


69 


3.129 


- 3.171 


0.598 - 


0.606 


0.192 - 


0.198 


0.8 


200.76 


0.22 


0.133 


- 0.139 


64 - 


65 


3.120 


- 3.162 


0.684 - 


0.693 


0.184 - 


0.190 


0.9 


210.97 


0.21 


0.127 


- 0.132 


61 - 


62 


3.111 


- 3.153 


0.769 - 


0.780 


0.177 - 


0.183 


1.0 


222.28 


0.20 


0.120 


- 0.125 


59 - 


59 


3.102 


- 3.144 


0.855 - 


0.866 


0.170 - 


0.176 



The masses obtained for both white dwarfs are un- 
usually high indicating that the star probably reached the 
asymptotic giant branch (AGB) in its evolution. The ini- 
tial binaries must have been very wide in order for this to 
happen. We notice that the masses calculated here differ sig- 
nificantly from those given in Section. 4, measured by fitting 
the line profiles to stellar atmosphere models. This discrep- 
ancy suggests that the redshift measurements may not be 
reliable, perhaps not surprising given the difficulty of sepa- 
rating the M star emission from the white dwarf absorption. 
Measurements of the white dwarf at UV wavelengths would 
be helpful. 

As mentioned in Section 3.3.1, the flux modulation seen 
in Fig. 7 for WD2009+622 indicates that the M dwarf com- 
panion is strongly irradiated and therefore the value mea- 
sured for Km is probably a lower limit for the true radial 
velocity semiamplitude. This implies that the M dwarf mass 
and inclination given in Table 8 are upper and lower limits 
respectively. To calculate how distorted the M dwarf is, we 
calculate how much the radial velocity semiamplitude of the 
companion changes as a function of its radius and how its 
mass and inclination are affected. We present the results in 
Table 9. The radius of the M dwarf is given in terms of a lin- 
ear filling fraction, /, defined as the ratio of the stellar radius 
measured from the centre of mass to the inner Lagrangian 
point. A value of / = 1 implies that the M dwarf fills its 
Roche lobe. For these calculations we have taken Mwd = 
0.61 and 0.64 Mq (the maximum and minimum values cal- 
culated above) and Kwd = —43.81 kms -1 . If the M star 
does not deviate too far from the main sequence we expect 
it to fill at least 0.4 of its Roche lobe which translates into a 
true radial velocity semiamplitude in the range 168 < Km < 
222 kms~\ and a mass between 0.120 < Mm < 0.165 Mq. 



3.3.3 The masses of the M dwarf companions 

An independent estimate of the masses of the M type com- 
panions can be done from their absolute infrared magni- 
tudes. Although the white dwarfs dominate the flux in the 
optical range in both systems, they produce only a minor 
fraction of the infrared luminosity. The distances of the 
systems were computed from the parameters of the white 
dwarfs given in Section 4. Since the flux contribution of the 
M dwarfs in the blue part of the spectra, used for the model 



atmosphere fits, is very small, we do not expect systematic 
effects caused by spectral contamination. 

J, H, and K magnitudes were retrieved from the 2MASS 
point source catalogue. We computed the white dwarf's con- 
tribution using the colour calibration of Bergeron et al. 
(1995) and subtracted it from the observed fluxes. Correc- 
tions are small for WD1042— 690 and do not exceed 25% 
(10%) for the J (K) band flux of WD2009+622. Finally, 
the M dwarf masses were computed from the calibration 
of Henry & McCarthy (1993). Results are listed in Ta- 
ble 10. Our error estimates include photometric errors and 
distance uncertainties resulting from the spectral analysis of 
the white dwarf. We adopted Napiwotzki, Green & Saffer's 
(1999) estimates for the external fit errors. 

For both systems, the value obtained for the mass of 
the M dwarf companion from its infrared magnitudes is in 
agreement with that obtained from the spectral analysis in 
section 3.3.2. In the case of WD2009+622, the revised mass 
after taking into account heating by the white dwarf is con- 
sidered. 



4 DISCUSSION 

Table 11 gives a list of all the detached white-dwarf binaries 
with known orbits known up to now. Most values have been 
taken from Ritter & Kolb (2003). The 6 systems discussed 
in this paper are also included. 

Liebert, Bergeron & Holberg (2005), Bergeron et al. 
(1992) and Bragaglia et al. (1995) obtain the temperatures 
and gravities for the systems studied in this paper. Their 
values, together with our determination for the masses of 
the white dwarf, are given in Table. 12. 

Fig. 9 presents the theoretical mass versus orbital pe- 
riod distribution for DDs from the population model as de- 
scribed in Nelemans et al. (2004). The four DDs discussed 
in this paper are plotted over the theoretical distribution 
and fall right in the expected range of values according to 
Nelemans et al. (2004). 

The direct progenitors of the double white dwarf bina- 
ries (giant plus white dwarf) could have had a wide variety 
of masses and periods, leading to inferred efficiencies of the 
CE that are poorly constrained (see Nelemans & Tout 2004, 
fig. 5). 

We calculated the possible progenitor systems of the 
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Table 10. Infrared properties and mass estimates for the M dwarf companions. 

dist (pc) Mj M H M K M/Mq 

WD1042-690 34.8±2.5 8.71±0.16 8.18±0.16 7.85±0.16 0.169±0.010 

WD2009+622 115±9 9.28±0.18 8.80±0.18 8.39±0.18 0.136±0.009 



Table 11. List of all the detached white dwarf binaries with known orbital periods (given in days). The type of binary is also given where 
WD = white dwarf; M = M dwarf; sdO/sdB = O/B subdwarf; ? = uncertain. * indicates periods measured in this paper. References for 
the orbital periods not measured in this paper arc (a) Bragaglia, Grcggio & Rcnzini 1990, (b) Koen, Orosz & Wade 1998, (c) Orosz & 
Wade 1999, (d) Maxted et al. 2000a, (e) Morales-Rueda et al. 2003a, (f) Maxtcd ct al. 2000b, (g) Marsh 1995, (h) Marsh et al. 1995, (i) 
Moran ct al. 1999, (j) Saffer, Livio & Yungclson 1998, (k) Napiwotzki ct al. 2002, (m) Maxted, Marsh & Moran 2002, (n) Holberg et al. 
1995, (o) Drechsel et al. 2001, (p) Kilkenny ct al. 1998, (q) Maxted et al. 1998, (r) Wood & Saffer 1999, (s) Orosz ct al. 1999 (t) Bruch & 
Diaz 1998, (u) Gizis 1998, (v) Wood, Harmcr & Lockley 1999 (w) Delfosse ct al. 1999 (x) Napiwotzki ct al. 2001, (y) Maxtcd ct al. 2000c, 
(z) Saffer, Liebert & Olszewski 1988, (aa) Hcbcr ct al. 2003, (ab) Karl ct al. 2003, (ac) O'Donoghue ct al. 2003, (ad) Maxted ct al. 2002, 
(ae) Hillwig et al. 2002, (af) Maxtcd ct al. 2004, (ag) Kawka ct al. 2000, (ah) Kawka ct al. 2002, (ai) Saffer ct al. 1993, (aj) O'Brien, 
Bond & Sion 2001, (ak) Sing et al. 2004, (al) Moran, Marsh k, Bragaglia 1997, (am) Morales-Rueda et al. 2003b, (an) Edclmann, Hcbcr 
& Napiwotzki 2002, (ao) Napiwotzki et al. 2004, (ap) Fuhrmeister & Schmitt 2003, (aq) O'Toole, Heber & Benjamin 2004, (ar) Robb 
& Greimel 1997, (as) Heber et al. 2004, (at) Raymond et al. 2003, (au) Gansicke et al. 2004, (av) Wood, Robinson & Zhang 1995, (aw) 
Bruch, Vaz & Diaz 2001, (ax) Pigulski & Michalska 2002, (ay) Shimansky, Borisov & Shimanskaya 2003, (az) Rauch & Werner 2003, 
(ba) Chen ct al. 1995, (bb) Green, Richstone & Schmidt 1978, (be) Landolt & Drilling 1986, (bd) Bell, Pollacco & Hilditch 1994, (be) 
Pollacco & Bell 1994, (bf) Lanning & Pesch 1981, (bg) Bleach ct al. 2002, (bh) Vcnncs & Thorstensen 1994. 
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Figure 9. Left panel: mass of the bright white dwarf for the four DDs discussed in this paper as measured by Bragaglia et al. (1995), 
Liebert et al. al. (2005) and Bergeron et al. (1992) (see Table 5) versus period distribution. In the grey scale we plot the mass to period 
distribution of DDs according to theory (for the model described in Nelcmans et al. (2004)). Right panel: mass ratio (only upper limits 
taken from Table 5) for the four DDs studied versus period distribution. The theoretical distribution according to the Nelemans et al. 
(2004) is also plotted. 



Table 12. Temperatures and gravities measured for the bright 
component of the system by fitting the hydrogen line profiles to 
stellar atmosphere models, (a) Bragaglia et al. (1995), (b) Liebert 
et al. (2005) and (c) Bergeron et al. (1992). V represents the 
V magnitude taken from the literature and My is the absolute 
magnitude. The masses given have been determined using cooling 
tracks by Althaus & Bcnvcnuto (1997). 
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two white dwarf plus M star binaries. If the high masses in- 
ferred from the gravitational redshifts were right, this means 
the direct progenitors of the white dwarfs must have been 
highly evolved giants. Using the equations in Hurley, Pols 
& Tout (2000) in the same way as described in Nelemans & 
Tout (2004) we calculated the possible progenitor systems. 
For WD1042— 690 we find possible progenitors with masses 
typically in the range 2 - 3.5 M Q , while for WD2009+622 
the progenitors typically have masses between 1.25 - 3 Mq. 
Because of the rather extreme mass ratios and the large 
radii of the giants, the CE is most likely caused by tidal in- 
teraction, rather than Roche-lobe overflow (Hut 1980) and 
the binaries will generally not be synchronised. We use the 
formalism derived for star - planet interactions by Soker 
(1996) to calculate the separation between the two stars at 
which the CE sets in. The required CE efficiencies are be- 
tween 0.1 and 1 for WD1042-690 and between 0.1 and 1.6 
for WD2009+622. Both systems can also be explained with 
the gamma-algorithm (Nelemans & Tout 2004), with values 
of 7 around 1.5. If, on the other hand, the lower white dwarf 
masses presented in Table 12 are right, the progenitor masses 
inferred are in the range 0.9 - 2.5 M for WD1042-690 and 



0.75 - 2 M for WD2009+622. The required CE efficiencies 
derived in this case lie in the same ranges as for the more 
massive white dwarf alternative presented above. 

Using the binary results presented in Tables 7 and 8, 
and the equations from Schreiber & Gansicke (2003), wc 
have calculated the white dwarf cooling age using models by 
Wood, Robinson & Zhang (1995), t coo i, the orbital period 
of the binary at the end of the CE phase, Pce, the time it 
will take for the binary to start mass transfer (and become a 
CV) assuming classical magnetic braking (CMB) and assum- 
ing reduced magnetic braking (RMB), t s( j, and the orbital 
period when mass transfer starts, P s d- These values are pre- 
sented in the top three rows of Table 13. The input masses 
used for the white dwarfs are those obtained assuming that 
they are carbon-oxygen core white dwarfs. The input masses 
for the companions are O.17M for WD1042-690 and 0.120 
< M M < 0.165 M for WD2009+622 (as calculated in Sec- 
tion 3.3.2). The value of Pce for both binaries is very close 
to their present orbital periods which indicates that they 
are very young PCEBs. These systems will not become CVs 
within a Hubble time (except perhaps WD1042-690 if RMB 
takes place), assuming To = 1.3 xlO 10 yrs (Ferreras, Mel- 
chiorri & Silk 2001), and when they do, their orbital periods 
will place them below the CV period gap. 

The evolutionary properties of WD1042— 690 and 
WD2009+622 are also calculated by using the lower white 
dwarf masses presented in Table 12 and the masses of the M 
dwarfs obtained in Section 3.3.3. These are presented in the 
bottom two rows of Table 13. The values for all parameters 
in this case are similar to those obtained for the larger white 
dwarf masses and the conclusions reached are equivalent. 



5 CONCLUSIONS 

We have obtained the orbital solution for four DD systems 
and two white dwarf - M dwarf binaries. We find that the 
white dwarf companions for the four DDs studied contribute 
between 10 and 20 per cent of the total luminosity and 
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Table 13. Evolutionary properties of WD1042-690 and 
WD2009+622. The top three rows present the results obtained 
assuming the white dwarf masses given in Table 8 (0.72 and 
0.59 Mq respectively). 1 and 2 indicate calculations obtained for 
WD2009+622 by assuming M M = 0.12 and 0.165 Mq respec- 
tively. The bottom two rows present the results obtained when 
the lower white dwarf masses given in Table 12 arc assumed in- 
stead. In this case the M dwarf masses calculated in Section 3.3.3 
are used. The values of i coo i and t s< j given in the table are in fact 
the logs of i coo i and t st j in years. 

WD t cool P CE (d) t sd P sd 

CMB RMB CMB RMB (d) 

1042-690 7.92 0.3381 0.3375 10.30 10.08 0.073 
2009+622 1 7.31 0.7412 0.7411 11.43 10.96 0.055 
2009+622 2 7.31 0.7412 0.7411 11.30 10.78 0.071 
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1042-690 7.65 0.3376 0.3372 10.39 10.09 0.073 
2009+622 7.20 0.7412 0.7410 11.43 10.87 0.061 

remain undetected. The masses and periods obtained for 
these systems agree with theoretical mass period distribu- 
tions (based on Nelemans et al. (2004)). 

In the case of the white dwarf - M dwarf binaries we have 
been able to measure the motion of both components. We 
find that there are signatures of strong irradiation of the sur- 
face of the M dwarf component in WD2009+622. The white 
dwarf masses calculated from their gravitational redshift are 
unusually high compared to those measured by fitting their 
hydrogen lines with stellar atmosphere models, and could be 
the result of the evolution of giant stars with masses between 
1.25 and 3.5 M Q that went through a CE phase as a result of 
tidal interaction. These two binaries are young PCEBs that 
will evolve into CVs, although not within a Hubble time, 
with orbital periods below the period gap. 
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